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ABSTRACT 

Context. The disk atmosphere is one of the fundamental elements of theoretical models of a protoplanetary disk. However, the direct 
observation of the warm gas (»100 K) at large radius of a disk (»10 AU) is challenging, because the line emission from warm gas 
in a disk is usually dominated by the emission from an inner disk. 

Aims. Our goal is to detect the warm gas in the disk atmosphere well beyond 10 AU from a central star in a nearby disk system of the 
Herbig Be star HD 100546. 

Methods. We measured the excitation temperature of the vibrational transition of CO at incremental radii of the disk from the central 
star up to 50 AU, using an adaptive optics system combined with the high-resolution infrared spectrograph CRIRES at the VLT. 
Results. The observation successfully resolved the line emission with 071 angular resolution, which is 10 AU at the distance of 
HD 100546. Population diagrams were constructed at each location of the disk, and compared with the models calculated taking into 
account the optical depth effect in LTE condition. The excitation temperature of CO is 400-500 K or higher at 50 AU away from the 
star, where the blackbody temperature in equilibrium with the stellar radiation drops as low as 90 K. This is unambiguous evidence 
of a warm disk atmosphere far away from the central star. 
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1. Introduction 

A warm disk atmosphere is a natural con sequence of a disk 
externally he a ted by stellar radiation (e.g. ICalvet et all I199U 
Aika wa et al.L |2002|) . Together with a somewhat cooler mid- 
plane and a disk inner-rim, it is one of the fundamental ele- 
ments of the current theoretical models on which our under- 
stand ing of a protoplanetary disk relies (e.g. iDullemond et al.L 
I2001I) . The emission from warm gas (»100 K) in a protoplan- 
etary disk h as been be s t observed in the vibrational transitions 
of CO (e.g. iGoto etail I2006L IBrittain et al.L I2007L iNaiita et al. , 
2007) and other molecules (e . g. IWe intraub. Kastner, & Barv, 
2000; ICarr. Tokunaga. & Nalitall2004l) . It is hard to estimate the 
physical properties of the disk atmosphere from this spatially 
unresolved line emission alone, however, because the emission 
from the warm gas is dominated by the s trong radiation from 
the hot inner disk closer to the central star (iBrittain et al.L [2007; 
van der Plas et al ., 2009). Sub-mm spectroscopy does probe the 
molecules in the outer disk (>100 AU) through the rotational 
transitions , but mostly cool gas of th e temperature 15-30 K (e.g.; 
iDent et alll2005blPanic et alll2010h . In order to measure the hot 
gas in the disk atmosphere at radii of 10 to 100 AU, one requires 
(1) high angular resolution to exclusively observe the faint emis- 
sion from the outer disk without being affected by the strong 



emission from the inner disk, (2) a proper spectroscopic probe 
that is sensitive to the warm disk atmosphere, many of which 
are in the infrared regime, and (3) the multiple transitions of 
that spectroscopic probe to measure the gas temperature quan- 
titatively. 

The goal of our observation is to measure the rotational ex- 
citation temperature of CO fundamental band at 50 AU away 
from a star, using an adaptive optics system to isolate the emis- 
sion of the outer disk from that of the inner rim. The near- 
est Herbig Be star to t he solar system, HD 100546 (B9.5Vne; 
iHouk & Cowlevl [1975^ is us ed as a te stbed, because of its close 
distance (d =103 pc; van den Anckerl H998) and high luminos- 
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away from the observers by 40 
with two_spiral arms seen at 200 



Panic & Hogerheiide, 2009) inclined 
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AU in the scattered light 



2005 ). Th e dust d isk is truncated at 10 AU 



* Based on data collected in the course of CRIRES program [084. C- 
0605] at the VLT on Cerro Paranal (Chile), which is operated by the 
European Southern Observatory (ESO). 



(IGradv et all 1200 111 

and inward (Bouwmanetal 1 12003b lAcke & van den Anckerl 
2006). The gas and the dust at the inner rim are warm, 200 K 
and 1000 K, respectively, under the direct irradiation o f the cen- 
tral star dBouwman et all I2003L IBrittain et all 120091) . The vi- 
brational trans i tion of CO and H2 have been spatially resolved 
(Britt ain et all I2009L Ivan der Plas et all 12009b ICarmona et ail 
2011), but the disk atmosphere has not been studied in depth, 
because the emission from the outer disk was not properly sep- 
arated from that of the inner rim. In this paper, we will follow 
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Table 1. Extraction apertures and excitation temperatures. 
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-48.7 (35.4) 
-4.9(±5.0) 


-2.5 (2.0) 
-22.1(17.7) 
-6.3(±5.0) 


-1.0(1.0) 
-8.9 (8.9) 
-6.0(±5.0) 


0.0(1.0) 
0.0 (8.9) 
0.0(±12.0) 
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<o) Location of the aperture center at the spectral extraction. The positions are given with respect to that of the central star in the unit of detector 
pixel. The numbers in parenthesis are the sizes of the apertures. Positions and sizes of the apertures in pixel unit are converted to the physical 

scale assuming the distance to HD 100546 being d =103 pc, and the pixel scale of CRIRES 86 mas. (c) Line center and full interval of the velocity 
range that the intensity of the line emission is integrated to calculate the line flux. Velocity offset is given with respect to the line center at the 

aperture location P4. <rf * The measurements are from the fittings of the population diagrams with LTE models presented in the red lines in Figure[2] 



up these studies, and measure the gas temperature in the disk 
atmosphere quantitatively. 

2. Observations 

The observation was carried out with CRIRES dKaufl et all 
2004) at the VLT on 29 and 30 April in 2010 in the open- 
time program 084.C -0605. The adaptive optics system MACAO 
(Bonnet et al., 2004) was used to feed nearly diffraction-limited 
images to the spectrograph using HD 100546 as a wave- 
front reference. The slit width was 072, providing spectra with 
a velocity resolution of 3 km s (R= 100,000). The spec- 
troscopy covers the wavelength from 4.588 pm to 5.004 pm 
[R(16)-R(0), P(l)-P(27)] continuously with six grating set- 
tings. The spectra were obtained with the observing tem- 
plate CRIRES spec jobs SpectroAstrometry. The slit was ori- 
ented eight position angles to increase the spatial coverage of 
the disk. We here used the subset of the data with the slit aligned 
to the major axis o f the disk projected onto the sky (P.A.=145°; 
lArdilaet al.ll2007l) . and anti-parallel to it (P.A =325°) to test the 
consistency of the detected line emissions, and the measurement 
of the excitation temperature of CO. The telescope was nodded 
by 10" every two exposures to subtract the thermal emission 
from the sky. The total integration time was 4 minutes for one 
grating setting. A spectroscopic standard star HR 6556 (A5 III) 
was observed at close airmass with HD 100546 with the same 
instrument settings. 

3. Data reduction and analysis 

The spectral images were pre-processed for subtraction of 
the dark-current images and flat-fielding, and coadded by 
cries_spec_jitter recipqj on the ESO gasgano platforrrQ. Because 
all line images of v =2-1 look qualitatively similar, those lines 
that were not blended with strong telluric absorption lines, or 
with the other transitions of CO lines, were registered and added 
to a single line image to increase the signal-to-noise ratio, and to 
check the spatial extent of the gas emission (Fig. |T}. Continuum 
emission was interpolated from both sides of the line emission, 
and subtracted. The gas emission of the disk is clearly resolved, 



1 CRIRES Pipeline User Manual VLT-MAN-ESO-19500-4406. 

2 http://www.eso.org/sci/data-processing/software/gasgano/. 



showing Keplerian rotation profiles consistently in both line im- 
ages obtained with slit P.A.=145° and 325° with the low-velocity 
wings extending up to 50 AU. The disk around HD 100546 is 
known to be brighter at the so uthwestern side tha n the northeast- 
ern side in the scattered light (lArdila et all 12007). The asymme- 
try of the surface brightness has been attributed to the preferen- 
tial forward scattering of the dust grains in the disk. The rotation 
of the disk detected in CO v =2-1 with the southeastern part ap- 
proaching toward us is consistent with this disk in clination and 
the counter-clockwise w inding of the spiral arms dGradv et all 
120051: lArdila et all 120071) . 

One-dimensional spectra were extracted from the spectral 
images, setting the extraction apertures at the incremental dis- 
tances from the central star (the positions and the sizes of the 
extraction apertures are summarized in Table [TJ. The telluric 
absorption lines were removed by dividing by the spectra of 
the spectroscopic standard star after correcting the small mis- 
matches in the wavelengths, the airmass, and the spectral reso- 
lutions. The wavelength calibration was performed to match the 
telluric absorption line s to the atm ospheric transmission model 
calculated by ATRAN jLordl[T992h . 

Because the transitions to the ground level v =1-0 were af- 
fected by the telluric absorption lines of the same transitions, 
v =2-1 lines were used to calculate the column densities of 
the rotational levels from 7=0 to J =26 in v =2 vibrational 
state to construct the population diagrams. The column density 
at each J level Nj was calculated for the optically thin limit 
NjA u [hv(L> = Wyfy, where u> is the solid angle of the emitting 
area, and f v is the flux density of the continuum emission to con- 
vert the equivalent width W v to the line flux. We were unable to 
perform the absolute flux calibration of the spectra to measure f v 
properly because of significant systematic errors from the vari- 
able PSF, variable slit transmission, and the variable pixel sam- 
pling of the PSF combined together. The factor f v /a> is irrelevant 
here, though, as long as it changes only slowly with the wave- 
length, because we are not primarily interested in the absolute 
column densities, but the relative distribution of Nj to measure 
the excitation temperature of CO. The line flux of the spectra 
obtained with adjacent grating settings were scaled so that the 
lines that were covered by both grating settings have the same 
line flux on average. The column densities divided by the statis- 
tical weights are shown in Fig.|2]after being shifted vertically by 
an arbitrarily chosen amount for the presentation. 
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Fig. 1. Top: Position-velocity diagrams of CO v=2-l from HD 100546. All vibrational transitions of CO v =2-1 in our observational 
coverage that do not overlap with the deep absorption lines of the terrestrial atmosphere are registered and added. The spectra were 
recorded with the slit P.A. = 145° (left) and 325° (right). The spatial axis is in abscissa with the position angles to which the slit was 
aligned indicated by arrows. The line image obtained with slit P.A.=325° is flipped horizontally so that the offset from the central 
star increases to the right in both of the line images. The continuum emission of the star was interpolated from both sides of the 
line emission and subtracted. The contours are drawn at every 10% of the peak brightness, except the lowest one which is at 5%. 
The apertures of the spectral extraction are marked by rectangles. Faint emission is still visible at 50 AU away at low velocity. The 
sizes and the locations of the extraction apertures are summarized in Table Q] Bottom: Spatial profile of the line emission, i.e., the 
position-velocity diagram shown in the upper panels are crushed along the velocity space. 



Obviously, the gas is warm (>400 K) out to the longest dis- 
tance (49 AU) covered by the observations. In order to measure 
the excitation temperature of CO quantitatively, model popu- 
lation diagrams were calculated based on synthetic line emis- 
sion spectra, and compared to the observations. The column 
densities in the upper and lower levels are given by N(„j) = 
N v= (i,2) g(uj)IQ(u,t) e~ E M /kT " in the LTE condition, where A^=(i,2) 
is the total CO column density at v = 1 or 2 state, g is the statisti- 
cal weight, Q is the partition function, £(„ /) is the energy of the 
upper or lower level with respect to the rotational ground level, 
and r ex is the excitation temperature. The line source function 
of a two-level system was given by S Y — N u A u i/(NiBi u - N u B u i), 
and the line emission spectrum was therefore calculated as I v = 
S v (l- e- T >) with the line opacity r v = g {N,B hl - N u B ul ) 0(v). 
The line profile 0(v) was assumed to be a Gaussian function. 
We performed preliminary runs to fit the observed population 



diagrams by the models with limited sampling of T ex and Nco, 

varying the line width of the Gaussian line profile from 2 km s _1 

to 30 km s _I . The fitting result in terms of the absolute deviation 

(Z |A^ bs - N?° del \) became better up to 6 km s _1 in most of the 
j 

spectra extracted from different disk locations. In several cases 
the fitting further improved up to 25 km s _! , but very slightly. For 
simplicity, we fixed the line width to 6 km s , and used the same 
line profile for the spectra extracted at all disk locations. The 
equivalent widths of the emission lines were measured with the 
model spectra and converted to the column densities in the same 
way as the observed spectra. The energy levels, the line cen- 
ter wavelengths, and Einstein A u i coefficients were taken from 
iGoorvitchl (119941) . A more detailed discussion on the calculation 
can be found elsewhere dGoto et al.Ll201 ll) . 
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Fig. 2. Population diagrams of CO v =2-1 based on the spectra extracted at the incremental distances from the central star. The 
upper and the lower panels are from the spectra recorded with the slit P.A.=145° and 325°. Open circles and filled triangles are for 
R- and P-branches, respectively. The column densities divided by the statistical weights are vertically shifted by arbitrarily amount 
for the presentation. The uncertainty of the individual column density is on the order of the dispersions of the data points. The red 
lines are the best-fit LTE model population diagrams assuming Gaussian line width of 6 km s _I . The blue lines are the same but the 
column density of CO is fixed to Nqo = 10 20 cirT 2 . 



The line emission is apparently optically thick with turnovers 
seen around Ej/k ~6500 K. The excitation temperature and the 
column density are heavily degenerated in these cases. The for- 
mal fitting error might underestimate the realistic uncertainty in 
the temperature measured. We produced contour maps of the ab- 
solute deviations between calculated and observed population 
diagrams on a grid of temperatures (500-2500 K) and column 
densities (10 19 -10 22 cuT 2 ) (Fig. IA.lllA.2l in the online edition). 
The total extents of the contour in temperature where the abso- 
lute deviation increases by 10% of the best fitting value were 
taken as the uncertainties of the excitation temperature of CO 
(Fig. [3] Table [JJ). We lack the measurements of the equivalent 
widths of high-/ levels at 49 AU because of the insufficient 
signal-to-noise ratio of the spectra (Fig. [2]). The upper limit of 
the excitation temperature is therefore less constrained at the ra- 
dius than at other locations, but the conclusion (T ex >400 K) is 
not affected. Fitting was also performed with the column densi- 
ties fixed to 10 20 cirT 2 to give an independent estimate of how 
the result is affected by this degeneracy. The resulting excitation 
temperatures do not differ by more than the extent of the errors 
estimated above except at -49 AU with P.A.=145°. 

4. Discussion 

The rotational excitation temperature of CO is between 400 K 
and 1100 K at all locations on the disk where the spectra were 
extracted. There might be a hint of the cooling of the gas down 
to 400-600 K with the distance from the central star, but this 
is not unambiguous. The measurements are consistent in the 
data with PA.= 145° and 325°, except at the edge of the disk 
at +49 AU, and the inner disk at +9 AU from the star. The latter 
might be picking up the emission from the hot inner rim ex- 
clusively. Although the observation was made with high angu- 
lar resolution, the line emission from different radii of the disk 
overlap up to 22 AU, because the emission from the inner rim 
is by far stronger than that from the disk behind (the spatial 
profiles of the line emission are shown in the lower panels in 
Fig. [TJ|. However, as is seen in the locations of the extraction 
apertures overlaid with the line images (Fig. [TJ, the spatial res- 
olution is sufficiently high that little blending with the emission 



from the inner disk is expected at 50 AU and beyond. The exci- 
tation temperature is still hotter than 400 K at 49 AU, where 
the blackbody equilibrium temperature with the stellar radia- 
tion [cr7± = LJ(16nd 2 ); cr, L„, d are the Stefan-Boltzmann 
constant, the luminosity of the star, and the distance from the 
star to the particle] drops as low as 90 K. Although the pres- 
ence of this warm ga s is exactly what the disk mo del predicts 
dKamp & Dullemondl 12001: Uonkheid etail 120041) . this is the 
first time that the temperature of the disk atmosphere is directly 
measured beyond the disk inner rim by a high angular resolu- 
tion observation that spatially resolved the disk. The range of 
the excitation temperature nicely overlaps with 300-800 K of 
the pure rotational lines of CO at high-/ (14< / <30) measured 
in the recent spectroscopy of the sam e object by Herschel/PACS 
(ISturm et all l2010h . iThi et al.1 (1201 lh also detected CFF /=5-4, 
6-5, 3-2 by the same instrument, and derived a rotational excita- 
tion temperature 323^ 2 | ° K. It is inferred that these observations 
are likely looking at the gas in a similar l ocation of the dis k as the 
present study. The model presented by IThi et alJ d201 lh . which 
takes into account the specific disk geometr y of HD 100546 con - 
strainedby the spectral energy distribution (Beni stv et al.L l2010). 
also shows a warm layer of the gas at 50 AU of the disk. 

A mechanism that likely vib rationally excites CO to v =2 
level is the ultraviolet pumping (IRrotkov et a because 
the v ibrational temperature (TVa, = 6600 ± 700 K; Ivan der PiasL 



2010 ) is much hotte r than the rotational excitation temperature 
(Brittain et al., 2009). The rotational levels are much faster to 
thermalize through the collisions with the molecular/atomic hy- 
drogen and therefore should represent the ambient gas tempera- 
ture better. It is in principle possible to account for the extended 
line emission not by the emission arising on the site, but by the 
scattered line emission, where CO is excited to v =2 level in the 
vicinity of the central star, and the line emission emerging from 
the hot inner disk is simply reflected into the line of sights to ob- 
servers by the dust grains at 50 AU. We discarded the scattered 
line emission, because the line images in Fig.[TJclearly show the 
Keplerian profiles, which indicate that the line center velocity 
at 50 AU is closer to the system velocity than that of the gas 
orbiting closer to the central star. 
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Fig. 3. Radial profile of the rotational excitation temperature of 
CO v =2-1 in the disk atmosphere of HD 100546 system. Data 
points are horizontally shifted slightly for clarity. The excitation 
temperatures are measured from the population diagrams pre- 
sented in Figure [2] The red squares are for an LTE model with 
fixed line width of 6 km s _1 . The line emissions are optically 
thick, therefore the column density and the excitation tempera- 
ture are highly degenerated. The error bars are given not by the 
fitting error in Figure [2] but by the total extent of the contour in 
temperature in the absolute deviation plots in Fig. IA.ll and lA.2l 
where the absolute deviation becomes worse by 10 % of the best 
fitting values. The blue triangles are from the fitting of the pop- 
ulation diagrams, but with the column density of CO fixed to 
N co = 10 20 cm" 2 . 



The observed exc itation temperature agrees with that of 
iJonkheid et al.1 (120071) . who calculated the chemistry and the 
temperature structure of a disk around a Herbig Ae star with- 
out the a priori assumption that the gas is physically coupled 
to the dust. The calculations were performed for the disk mass 
from 10 _1 to 1O~ 4 M to simulate the dissipation of a disk, and 
with and without the dust settling, i.e., the mass ratio of the dust 
and the gas in the upper layer gradually reduced from 0.01 to 
10~ 5 . The gas temperature in the upper layer is clearly higher in 
cases without dust settling, because the gas is primarily heated 
by the photoelectric effect of the dust grains. The gas is hotter 
than 300 K in the disk atmosphere at the distance of 50 AU from 
the disk surface to the significant depth toward the mid-plane 
for all the disk models calculated without dust settling, which is 
comfortably matched with the excitation temperature we found; 
while the gas temperature with dust settling is everywhere less 
than 100 K. That the disk atmosphere is well mixed with the dust 
grain s is also cons i stent with the classification of the source as 
la by lMeeus et al.l (120011) with spatially resolved (jHa bart et al., 
2006), rich emission features of PAHs (lAcke et all l2010h . and 
the thermal continu um emission extended up to 1'.'4, or 145 AU, 
in the mid-infrared (Mu lders et al.Ll20Tlh . 
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Appendix A: Measurement of rotational excitation 
temperature 

Here we provide the contour maps of the absolute deviations 
between calculated and observed population diagrams that we 
used to estimate the uncertainty of the excitation temperature 
of CO. The extents of the contour in the temperature where the 
absolute deviation increases by 10% of the best fitting value were 
taken as the error of T ex . 
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Fig. A.l. Contour maps of the absolute deviations between the calculated and observed population diagrams. The observed popu- 
lation diagrams are based on the data obtained with the slit P.A.=145°. T ex and Nco that give the minimum absolute deviations are 
marked by crosses. The contours are drawn at every 10% increase of the absolute deviation from its minimum value. 
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Fig. A.2. Same as Fig. IA.ll but for the data obtained with the slit P.A.=325°. 



